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The method involves preliminary triangulation for 
modelling of boundaries, and division of the 
region into blocks with up to six faces 
corresponding to discontinuities or principal 
faults. Each face is reticulated by interpolation 
from its edges to form an intermediate grid, 
whose nodes are projected from a defined 
direction on to the triangulated surface. Each 
block is stocked with a regular distribution of 
internal nodes by interpolation from those of the 
grids formed on the various faces. The 
positioning of nodes of the elementary cells is re- 
arranged by relaxation in accordance with 
predefined criteria. 
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Description 

[0001] La presente invention concerne une methode pour g§n6rer un maillage 3D dans le but de modeliser un corps 
queiconque dont la geometrie est connue et imposee a priori tel que par exemple une structure geologique en tenant 

5 compte fidelement de sa geometrie dans toute sa complexite naturelle. Le maillage realise par la methode s'accorde 
bien dans le cas d'une structure geologique, avec la position des interfaces entre les couches, avec les failles meme 
de forme complexe, les biseaux et autres irregularites de structure, en les representant avec leur conformation reelle. 
[0002] La methode selon P invention trouve des applications notamment pour la modelisation en geomecanique, pour 
la caracterisation des reservoirs : couches de sediments, stratigraphie, et aussi pour la simulation des ecoulements 

10 dans les bassins et reservoirs. 

I) ETAT DE LA TECHNIQUE 

[0003] La modelisation d'une structure geologique est realisee a partir cfun ensemble de donnees disponibles ac- 
15 quises soit par des mesures, so it par des analyses et des interpretations de donnees sismiques ou des diagraphies 
realisees dans des puits. Ces donnees se presentent sous ta tonne d'un nuage de points entre lesquels, au mieux, it 
convient d'installer un maillage. La formation d'un maillage en 3D joue un r6le crucial pour faire des simulations nu- 
meriques dans differentes applications geologiques. Beaucoup de recherches ont ete menses sur le sujet principale- 
ment en ingenierie de reservoir en raison de ses applications dans le domaine petrolier. 
20 [0004] On connait de nombreuses techniques de maillage 3D pour des simulations numeriques par elements ou 
differences finis et des ensembles de logiciels de simulation crSes pour leur mise en oeuvre qui sont adaptes a des 
applications et respectent des contraintes determinees. 

[0005] II est connu pour certaines applications d'utiliser des maillages de type tetraedrique non structure. II est connu 
aussi de faire du maillage hexaedrique de peau ou de surface et de repartir des points dans le maillage soit par 

25 interpolation transfinie, soit encore par resolution d'equations aux derivees partielles. Les surfaces sont generalement 
creees a partir de fonctions parametriques et les points du maillage sont obtenus simplement par echantillonnage de 
ces fonctions. Une telle methode ne convient pas generalement pour des applications dans ies geosciences ou les 
volumes a modeliser ont des formes le plus souvent tres complexes, difficiles a traduire par des fonctions parametriques 
compte-tenu de la precision qui est recherchee. Aussi prefere t-on approcher les formes complexes par des elements 

30 de surface trianguies. 

[0006] Differentes techniques de maillage appliquees a la modelisation de reservoir, sont decrites par exemple dans : 

Georges P.L. : "Generation Automatique de Maillages", Editions Masson, 1991 ; 

35 - Gugrillot D.R. et al: An ineractive 3D Mesh Builder for Fluid Flow Reservoir Simulation; in SPE - Computer Appli- 
cations, Dec. 1993; 

Aziz K.; Reservoir Simulation Grids : Opportunities and problems. SPE-JPT, July 1993; ou 

40 - Heineman Z.E. et al : Gridding Techniques in Reservoir Simulation. Intl Forum on Reservoir Simulation. Atpbach 
Austria, Sept. 12-16, 1988. 

[0007] Cependant, la plupart des techniques de maillage connues, et de la les simulations qui en resultent, ne con- 
viennent bien que pour modeliser des structures geologiques relativement simples. Les irregularites geologiques cou- 

45 rantes dans les zones-reservoir modeiisees, telles que les couches ayant de grandes variations cTepaisseur, les failles 
ou les biseaux, sont rarement pris en compte ou alors seulement sous une forme grandement simplifiee. Le document 
RENARD P ET AL: 'Three-dimensional geometric modeling of a faulted demain: the Soultz Horst example (Alsace, 
France)' COMPUTERS & GEOSCIENCES, NOV. 1994, UK, vol. 20, no. 9, pages 1379-1390, prend en compte les 
irregularites geologiques en prevoyant un traitement sp6cifique des volumes des couches presentant des irregularites 

so par construction d'une surface trianguiee qui passe au mieux par les donn6es existantes et dont la g6ometrie s'accorde 
avec la configuration geologique. 

[0008] Dans les modeles mailtes obtenus dans le cadre des geosciences, les failles apparaissent verticalement; au 
mieux, elles sont inclinees et planes. Or il est important de rendre fidelement, la geometrie en 3D des sediments car 
lis possedent des caracteristiques physiques propres, et aussi la g6ometrie des failles et autres discontinuites meca- 
55 niques. 

[0009] La methode selon 1'invention permet de gen6rer un maillage 3D respectant la g6om6trie d'une zone d'un 
milieu complexe dont la geometrie est connue et imposee a priori, tel qu'une structure geologique, dans le but de 
r6aliser un modfele representatif de ce milieu. La structure est detinie a partir de donnees disponibles acquises soit 
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par des mesures, soit par des analyses et des interpretations de donnees sismiques ou des diagraphies r6alisees 
dans des puits. Elle comporte une moderation prealable des limites (bords de la zone, discontinuites ou failles par 
exemple) par des surfaces triangulees, enveloppant au mieux ces limites, et un fractionnement de la zone d'elude 
suivant en blocs ayant au plus six faces, dont les faces correspondent a des discontinuites ou des failles principales 
de la zone prealablement triangulees. 
[0010] Elle est caract6ris6e en ce que'elle comporte : 

la subdivision de chaque bloc en un ensemble de cellules elementaires necessaires pour modeliser la zone par 
elements finis, par une reticulation en 3D de chacun des diffeYents blocs obtenus, subdivision qui comprend : 

a) le maillage de chacune des faces du bloc par interpolation a parti r de ses bords, de facon a cr§er une grille 
intermediaire, suivi d'une projection sur la surface triangulee des noeuds de cette grille intermSdiaire suivant 
une direction definie (perpendiculairement a elle-meme par exemple); 

is b) le peuplement en noeuds de grille de chaque bloc dans son 6paisseur, par une methode cfinterpolation a 

parti r des noeuds des differentes grilles formees sur les differentes faces triangulees de chaque bloc, et 

c) un rearrangement dans le positionnement des noeuds des cellules 6l6mentaires cr6ees, obtenu (de fagon 
iterative, g6ne>alement) par relaxation dans le positionnement des noeuds, en fonction de criteres predefinis 
20 tels que la regularite, le choix de la taille relative des differentes cellules et le mode de progression de leure 

tallies, la direction de projection etc. 

[0011] Chaque bloc etant ainsi reticule en 3D, on peut proceder a un rSassemblage ou reempilement de tous les 
blocs. 

25 [0012] On peut proceder successivement a une premiere operation de relaxation dans le positionnement des noeuds 
sur chacune des differentes faces triangulees, le peuplement en noeuds internes £tant effectug a parti r des grilles 
relax ees, et ensuite a une deuxieme operation de relaxation dans le positionnement de ces noeuds internes. 
[0013] De preference, la premiere operation de relaxation est effect uee par des phases iteratives de relaxation- 
projection. 

30 [0014] La mSthode selon Pinvention permet d'obtenir un modele qui respecte tres precisement la g6om6trie souvent 
complexe de structures geologiques telles que des champs petroliferes, ce qui a pour effet que les simulations qui 
peuvent etre faites a partir de ce modele rendent plus fidelement compte des ph6nomenes. La melhode utilised peimet 
en outre une economietres importante de temps de calcul pour realiser les operations de maillage 3D habituellement 
tres longues avec les methodes existantes. 

0015] D'autres caracteristiques et avantages de la m€thode selon invention, apparattront a la lecture de la des- 
cription ci-apres de modes de realisation d6crits atrtre d'exemples non limitatifs, en se r6f6rant aux dessins annexes ou : 



35 



40 



45 



les Fig. 1 et 2 montrent des exemples de relaxation-projection sur un maillage surfacique 4X4; 
la Fig. 3 montre un cube defonme avec des faces triangulees; 
la Fig.4 montre I'enveloppe exteme d'un cube maill6; 

tes Fig.5 et 6 montrent les surfaces internes d'un cube maill6 respectivement avant et apres relaxation; 

la Fig. 7 montre des courbes relevees sur cinq coupes d'une structure geologique de tres grandes dimensions, sur 
laquelle la m&thode selon Pinvention a 6t6 test&e, qui constituent une partie des donnees initiates dont on dispose, 
sur cette structure; 

so - la Fig. 8 montre d'autres donn6es initiales sur la meme structure, en I'occurrence, des cartes cfisobathes et les 
emplacements de 4 puits; 

la Fig. 9 montre une moderation surfacique ( base de surfaces triangulees) de la structure geologique 6tudiee; 
55 - (a Fig. 10 montre un exemple de decoupage de la structure en blocs; 

la Fig. 11 montre le resultat cfun maillage bloc par bloc de I'enveloppe exteme de la structure geologique; 
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la Fig. 12 montre le r6sultat (fun maillage bloc par bloc des surfaces internes de la structure geologique; 

la Fig. 13 iilustre a titre cfexemple comment on prend en compte un biseau dans le maillage de la structure geo- 
logique; 

5 

la Fig.1 4 iilustre a titre d'exempte comment on prend en compte une faille dans le maillage de la structure geolo- 
gique; 

la Fig. 15 iilustre la methode cTinterpolation 2D de Coons; 

10 

la Fig.1 6 montre les elements pris en compte dans la methode d 1 interpolation 3D de Coons; et 

la Fig.1 7 iilustre le mode de relaxation-projection des noeuds de mailie. 

is [0016] La structure geologique de grandes dimensions illustree par les Fig.7-10. correspond a un dome anticlinal 
comportant six interfaces entre des sediments. 

[0017] La realisation d'un maillage en 3D par elements finis d'une zone a partir d'un ensemble de donnees inrtiales 
connues pour la zone : interfaces entre diff6rentes couches ou horizons, failles principales, carte en courbes de niveau 
avec les emplacements de puits passe par une succession d'etapes qui seront decrites ci-apres. 

20 

I) Moderation des limites ou bords de la zone etudiee et de ses principales discontinues. 

[0018] On realise une modelisation des surfaces : bords, interfaces, failles en utilisant une logiciel specialise connu 
tel que GOCAD (marque deposee) sous la forme d'une juxtaposition continue d'elements de surface triangulaires. 

25 Dans un premier temps on forme une ebauche qui s'adapte au mieux aux donnees initiales. De maniere iterative, on 
deplace les sommets de la surface en utilisant un logiciel d'interpolation et on densifie localement les elements de 
surface triangulaires, de facon a plaquer au mieux la surface sur le semis de points (donnees d'entree). Les elements 
de surface ainsi crees ne se raccordent generalement pas comme ils le devraient : les interfaces ne sont pas limitees 
aux intersections et au bord du domaine d'etude et elles peuvent se croiser en debordant. On calcule alors les lignes 

30 d" interception entre surfaces et on limite celles-ci aux courbes d'intersection. La Fig 7 montre le resultat final de cette 
modelisation des limites de la structure et Ton peut remarquer que la plupart des triangles des surfaces triangulees du 
biseau sont localises dans les regions de forte courbure. 

[0019] Acestade, tl n'est pas encore possible deproduireun maillage 3D tenant compte des contraintes geologiques 
et numeYiques. Pour y parvenir, on procede a un fractionnement de la zone. 

35 

II) Fractionnement de la zone etudiee 

[0020] On decoupe d'abord la zone en macro-blocs ayant au plus six faces. II s'agit generalement de blocs de forme 
hexaedrique (a six faces), mais le n ombre de faces peut etre inferieur a six et descendre jusqu'a deux. Dans ce cas, 

40 pour les besoins de la methode, le bloc est quand meme consideYere comme un block hexaedrique mais ayant des 
faces degenerees (reduites a un point ou a une courbe). Le mode de decoupage est generalement suggere par la 
geometrie des horizons, des principales failles et leur interaction. Pour constituer ces macro-blocs 3D, on a besoin 
parfots d'introduire des frontieres artificielles. C'est necessaire par exemple quand un bloc a plus de six faces naturelles 
et ne satisfait done pas a la regie chotste d'etre au plus hexahedrique. Dans un tel cas, on reunit deux faces en une 

45 face complexe pour satisfaire le critere, ou bien on introduit des surfaces de coupure artificielles. 

[0021] Un exemple de decoupage de la structure geologique etudiee en macro-blocs est montre a la Fig. 1 0, ou Ton 
a extrapolS la gaille jusqu'a la frontiere de la structure et verticalement projete I'intersection du biseau sur les autres 
surfaces. 

so III) Maillage 

[0022] On realise une reticulation ou maillage en 3D s§parement de chacun de ces macro-blocs de facon a le sub- 
diviser en blocs 6l6mentaires qui vont former les cellules pour la modelisation par elements finis recherch6es. Cette 
operation est realisee en deux temps. 

55 

1 ) Maillage des faces 

[0023] a) on cree tout tfabord une grille sur chacune des cinq ou six faces de chaque macro-bloc. Pour cela, on 
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forme une grille intermediaire par interpolation a parti rdes quatre bords de la face, en utilisant la methode cfinterpolation 
transfinie de Coons dont le principe est rappele ci-apres. 

[0024] Etant donne quatre courbes parametriques f1(u), f2(u), g1(v), g2(v) (0 2 u, v 2 1) detinissant quatre bords 
jointifs (Fig. 15), on introduit une formule pour calculer la surface la plus tendue qui passe par les quatre bords (qui 
5 interpole ces bords) : 



5(^iO = (l-w)ffl(tO + i^2(v) + (l-v)/l(u) + v/2(M)-[(l-«)(l-v)P(0.0) 
10 +0 -i0vP(0.1) + «(1 -v>P(1.0) + wvP(l.l)] k 

[0025] Les sommets P(i, j) sont les quatre coins de la surface. 

[0026] On peut alors obtenir un maillage regulier NxM de cette surface par simple echantillonnage en prenant pour 
is sommets les points S(i, j) correspondant aux ui = i/(N-1 ), vj =j/(M-1), i=0 a N-1 , j=0 a M-1 . La formule 1 restreinte au 
points du maillage devient : 

20 5(f, y J = ( 1 - — i-j^Ky ) + ^ 2(7) + ( 1 - -j^IJ )/l (0 + -35fe /2 (/) 

"[(I " AM W " M^T + (1 - jfa PlO, I ) (2) 

25 

[0027] Dans le cas illustre, les courbes f 1 , f2, g1 , g2, sont obtenues a parti r du logictel GOCAD precite et sont done 
des lignes bris6es. Pour obtenir un maillage regulier MxN, on echantillonne regulierement (en abscisse curvtltgne ou 
longueur de segment) les courbes f 1 et f2 en N points, et les courbes g1 , g2 en M points, puis on applique la formule 
30 (2). On peut tout ausst bien obtenir un maillage a progression geom&trique sur chaque direction, simplement en im- 
posant cette contrainte au niveau du reechantillonnage des courbes. 

[0028] Ce maillage intermediaire sauf cas exceptionnel ne convient pas car il n'y a pas de raison qu'il coincide a ce 
stade avec la surface qui a prealablement ete triangulee par le logicie! de surfacage. L'ajustement necessaire consiste 
a projeter les noeuds de la grille intermediaire sur la surface triangulee orthogonalement a celle-ci, ou suivant cf autre s 

35 directions , de facon a obtenir une nouvelle grille avec des noeuds places sur la surface triangulee. Quand les quatre 
courbes frontieres sont regulierement echantiIlonn6es, la methode d'interpolation de Coons presente I'avantage de 
produire une grille reguliere en ce sens que sur chaque courbe isoparametrique tous les segments echantillonnes ont 
la m§me longueur. Plus g6neYalement la methode de Coons conserve les rapports de taille entre elements cons6cutifs 
dans une direction donnees. Ceci permet d f avoir par exemple des maillages reguliers ou respectant certain es pro- 

40 gressions dans la taille des mailles. Cependant la grille obtenue par cette methode connue perd sa regularity (ou le 
critere recherche) une fois qu'elle est projetee sur la surface triangulee. 

2) Rearrangement des noeuds de face 

45 [0029] Pour retro uver le critere desire dans la distribution des noeuds de la grille projetee, une regularity satisfaisante 
par exemple, il faut deplacer les noeuds du maillage sur la surface triangulee. Pource faire, on precede a une relaxation- 
projection de la grille par ajustements successifs comme indique ci-apres (dans le cas ou le critere est la r6gularite), 
en relation avec les Fig. 1 , 2. 

[0030] Etant donne un bloc mailie, le but cette relaxation est de rendre son maillage le plus regulier possible au sens 
50 ou n'importe quelle courbe isoparametrique du maillage soit regulierement echantillonn6e. Pour ce faire. on deplace 
de facon iterative tout point du maillage de maniere a le positionner au mieux a egale distance de ses deux voisins 
sur chacune des trois directions (u,v,w) de voisinage. Soit P le point qu'on cherche a deplacer. Pour chaque direction 
i, (i = u, v,w) (Fig.1 7), on trouve le point Pi qui se situe a egale distance des deux voisins de P sur la courbe de direction 
i. (i = u, v,w). Le point P est alors remplace* par le barycentre des trois points Pi (i = u, v, w). 
55 [0031 ] Apres avoir deplace tous les points du maillage de la sorte, ceux-ci ne sont plus sur la surface triangulee. On 
les projette alors a nouveau sur la surface triangulee et on repete le processus de relaxation -projection jusqu' a atteindre 
la regularity souhaitee. 

[0032] Ceci est illustre sur I'exemple des Fig. 1 , 2, correspondant a une face du cube detorme de la Fig. 3 qui est 
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un exemple synthetique permettant de bien Hlustrer les differentes etapes de la methode. Cette surface a ete obtenue 
par simple deformation d'une face carree en pnSservant ses bords ( le bord de la surface reste un carre). Sur ces 
figures, on peut voirtrois maillages 4X4. Le premier est un maillage plat issu de ("interpolation de Coons appliquee 
aux quatre bords (formant un carre). Le deuxieme est le maillage fin dont les noeuds sont obtenus par projection de 
s ceux du premier maillage. On constate que ce maillage a perdu de sa regularite. Le troisieme est le maillage regulier 
trace en gras et dont les noeuds sont situ6s sur la surface. II a ete obtenu a partir du deuxieme maillage apres une 
serie de projections et relaxations. La Fig. 4 montre le resultat de ces operations appliquees a I'ensemble des faces 
du cube deform e. 

10 3) Peuplement en noeuds de grille de chaque bloc dans son epaisseur : 

[0033] Pour prolonger en profondeur les grilles formees sur les faces des blocs (Fig.5), on utilise une generalisation 
au volume de la methode rfinterpolation transfinie precedente. 

[0034] Etant donne 12 courbes parametriques f1 (u), f2(u ), f3(u), f4(u), g1 (v), g2(v), g3(v), g4(v), h1(w), h2(w), h3 
is (w) f h4(w) (0 2 u, v, w 2 1) definissant les 12 bords d'un bloc tel que montre" dans la Fig. 16, une generalisation de la 
formule de Coons fournit une parametrisation volumique a I'interieur du bloc. La formule de Coons dans le cas volu- 
mique est la suivante : 



20 



V(u n v, vv) = (I - v»)( I - +(l - v)vv/2(w)+ vwf3(u) + v(| - w) f4(u) 

~(i - u )(l - w)gl{ v) +(!-«) wg2{v) + uwg3(v) + it( I - w)g4(v) 
■K l - «X I - v)h K*') + H - u)vh2{\v) + uvhX w) - «( I - v)h4(v) 

25 + C(«.l\H\) 



30 



35 



(3) 



avec 



C( u,v , w)=-2[(1 -u)(1 -v)(1 - w) P(0,0,0)+(1 - u)(1 -v) wP(0,0,1 ) 
+(1-Lf)v(1-w)P(0,1 l 0)+(1-u)vwp(0,1,1)+u(1-v)(1-w)P(1 l 0,0) 
-v) ivP(1 ,0,1 )+uv(1 -iv)P(1 ,1 ,0)+uv ivP(1 .1.1)] 



[0035] Les sommets P(i, j, k) (i, j, k= 0 ou 1 sur la Fig 1 6) sont les 8 coins du bloc. 

[0036] En combinant les equations (1) et (3) on aboutit a une expression de la formule de Coons dans le cas volu- 
mique en fonction des 6 surfaces de bord du bloc, S1(u, w), S2(u,w), S3(v,w), S4(v,w), S5(u,v), S6(u,v). Mais ici, les 
40 surfaces sont aussi supposees etre obtenues par interpolation de Coons a partir de leurs courbes de bord. La formule 
est la suivante : 



45 V(u y v, w) = -[(1 - v)S\(u y \v) + vSliu. vv) + ( 1 ~ u)S3(v,w) + m54(v. w) 

* 

+0 -w)S5(u,v) + vi/S6(ii % v))+-C(k,v.iv) 

4 

so 

[0037] Un echantillonnage quelconque de V(u,v,w) fournit alors un maillage EF du bloc. 

[0038] Cette interpolation en profondeur donnerait une distribution reguliere de noeuds si elle etait appliquee a partir 
des noeuds des grilles formes sur les differentes faces du bloc. Ce n'est pas le cas et la distribution obtenue est 
gen6ralement irr6guliere car la methode cTinterpolation en profondeur est appliquee aux noeuds reportes par relaxation 
55 sur les surfaces triangul6es, comme d6crit dans I'etape 2 precedente. 
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4) Rearrangement des noeuds de bloc 

[0039] Aussi, pour retrouver une regularite suffisante (ou respecter un autre critere choisi), la methode selon Pinven- 
tion comporte 6galement une relaxation de la grille 3D de la maniere indiquSe pr6c6demment. Les Fig.5, 6 montrent 
5 le maillage interne du cube d6forme avant et apres relaxation. Le maillage de la Fig. 6, on le note, possede une regularite 
bien meilleure que !e maillage de la Fig.5. 

[0040] Quand chacun des blocs precedents de la structure etudiee a ete convenablement maille en volume, on 
reassemble tous ces macro-blocs apres leur reticulation 3D. 

[0041] La methode qui vient d'§tre decrite peut etre utilisee pour mailler des h6t6rog§n6ites diverses telles qu'un 
10 chenal ou encore pour mailler radialement un puits devie. 

[0042] On a decrit des modes de realisation de la methode appliquee au maillage cfune structure geologique. On 
ne sortirait pas du cadre de I'invention, cependant, en utilisant la methode pour le maillage de tout corps ou structure 
dont la geomStrie est conn ue et imposee a priori. 

15 

Revendications 

1. Methode pour g6n6rer un maillage 3D dans une structure dont la g6om6trie est connue et imposed a priori telle 
que par exemple une structure geologique definie a partir de donnSes disponibles acquises par des mesures, par 

20 des analyses et/ou des interpretations de donnees sismiques ou des diagraphies realisees dans des puits au 

travers de la structure, dans ie but de realiser un mo dele representatif de cette structure, comportant une mod6- 
lisation prealabfe des limites par des surfaces triangu lees, enveloppant au mieux ces limites, et un fractionnement 
de la zone en blocs ayant au plus six faces, dont les faces correspondent a des discontinuites ou des failles 
principles de ia zone prealablement triangulees, la methode etant caracterisee en ce qu'elle comporte : 

25 

la subdivision de chaque bloc en un ensemble de cellules eiementaires necess aires pour modeliser la zone 
par elements finis, par une reticulation du maillage en 3D de chacun des differents blocs obtenus, subdivision 
qui comprend : 

30 a) le maillage de chacune des faces du bloc par interpolation a partir de ses bords, de facon a creer une 

grille intermedial re, suivi cTune projection sur la surface triangulee des noeuds de cette grille intermediate 
suivant une direction definie, de facon a generer des grilles de face; 

b) le peuplement de chaque bloc en noeuds internes par interpolation a partir des noeuds des differentes 
35 grilles formees sur les differentes faces triangulees de chaque bloc, de facon a obtenir une distribution 

reguliere des noeuds dans le volume de chaque bloc; et 

c) un rearrangement dans le position nement des noeuds des cellules eiementaires creees, obtenu par 
relaxation dans le positionnement des noeuds. en fonction de criteres predefinis. 

40 

2. M6thode selon la revendication 1 , caracterisee en ce que I'on procede success ivement a une premiere operation 
de relaxation dans le positionnement des noeuds sur chacune des differentes faces triangulees, le peuplement 
en noeuds internes etant effectue a partir des grilles relaxees, et ensuite une deuxifeme operation de relaxation 
dans le positionnement de ces noeuds internes. 

45 

3. Methode selon la revendication precedente, caracterisee en ceque Ton realise la premiere operation de relaxation 
par des phases iteratives de relaxation-projection. 

4. Methode selon Tune des revendications 1 ou 2 ou 3, caracterisee en ce Ton choisit comme direction de projection, 
so une direction perpend iculat re a la grille intermedial re. 

5. Methode selon Tune des revendications precedentes, caracterisee en ce que Ton utilise un mode d'interpolation 
de type transfini. 

55 6. Methode selon Tune des revendications precedentes, caracterisee en ce que Ton rearrange le positionnement 
des noeuds des cellules eiementaires creees de facon a obtenir des maillages reguliers. 

7. Methode selon Tune des revendications 1 a 5, caracterisee en ce que Ton rearrange le positionnement des noeuds 
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des cellules elementaires creees de facon a obtenir un mode determine de progression de la taille des mailles. 

8. Methode selon Fune des revendications prec6dentes, caracterisee en ce que Ton proofed e a un ^assemblage 
de tous les blocs. 



Patentanspriiche 

1 . Verfahren zum Erzeugen eines 3D-Netzwerkes in einer Struktur, deren Geometrie bekannt und a priori vorge- 
10 schrieben ist, wie zum Beispiel eine geologische Struktur, die aus verfugbaren Daten definiert ist, welche durch 

Messungen, durch Analysen und/oder Interpretationen seismischer Daten oder Diagraphien aufgenommen sind, 
welche in Bohiiochern entlang der Struktur realisiert sind, mit dem Ziel, ein fur diese Struktur representatives 
Modell zu verwirklichen, umfassend eine vorherige Modellierung der Grenzen durch aus Dreiecken zusammen- 
gesetzte Oberflachen, welche diese Grenzen am besten umhullen, und eine Fraktionierung der Zone in Blocke 
is mit hochstens sechs Flachen, deren Flachen Hauptdiskontinuitaten oder Hauptverwerfungen der zuvor aus Drei- 

ecken zusammengesetzten Zone entsprechen, wobei das Verfahren daclurch gekennzefchnet ist, dass es um- 
fasst: 

die Unterteilung jedes Blockes in eine Zusammenstellung von Elementarzellen, die notwendig sind, urn die 
20 Zone durch Endelemente durch eine Vernetzung oder 3D-Netzwerk von jedem der verschiedenen erhaltenen 

Blocke zu modellieren, welche Unterteilung umfasst: 

a) Das Netzwerk von jeder der Flachen des Blocks durch Interpolation ausgehend von seinen Grenzen 
derart, dass ein Zwischengrtter erzeugt wird, gefolgt von einer Projektion der Knoten dieses Zwischengit- 

25 ters auf die aus Dreiecken zusammengesetzten Fiache gemaB einer definierten Richtung, derart, dass 

Flachengitter erzeugt werden; 

b) die Population jedes Blockes mit inneren Knoten durch Interpolation aus Knoten von verschiedenen 
gebildeten Gittern auf die verschiedenen aus Dreiecken zusammengesetzten Flachen jedes Blocks, der- 

30 art, dass eine regelmassige Verteilung der Knoten in dem Volumen jedes Blockes erhalten wird; und 

c) Einer Wiederanordnung in derPositionierung der Knoten dererzeugten Elementarzellen, erhalten durch 
Entspannung in der Posrtionierung der Knoten in Abhangigkeit vorgegebener Kriterien. 

35 2. Verfahren nach Anspruch 1 , dadurch gekennzeichnet, dass man aufeinander folgend zu einem ersten Entspan- 
nungsvorgang in der Positionierung der Knoten auf jeder der verschiedenen aus Dreiecken zusammengesetzten 
. Flachen ubergeht, wobei die Population mit inneren Knoten ausgehend von den entspannten Gittern durchgefuhrt 
wird und man anschlieRend zu einem zweiten Entspannungsvorgang in der Positionierung dieser inneren Knoten 
ubergeht. 

40 

3. Verfahren nach dem vorhergehenden Anspruch, dadurch gekennzeichnet, dass man den ersten Entspannungs- 
vorgang durch iterative Relaxations -Projektions-Phasen durchfuhrt. 

4. Verfahren nach einem der Anspruche 1 oder 2 oder 3, dadurch gekennzeichnet, dass man ats Projektions- 
45 Richtung eine Richtung wahlt, die senkrecht zum Zwischengitter liegt. 

5. Verfahren nach einem der vorhergehenden Anspruche, dadurch gekennzeichnet, dass man einen Interpolati- 
onsmodus vom transfiniten Typ verwendet. 

50 6. Verfahren nach einem der vorhergehenden Anspruche, dadurch gekennzeichnet, dass man die Positionierung 
der Knoten der erzeugten Elementarzellen derart wiederanordnet, dass regelmaBige Netzwerke erhalten werden. 

7. Verfahren nach einem der Anspruche 1 bis 5, dadurch gekennzeichnet, dass man die Posrtionierung der Knoten 
der erzeugten Elementarzellen derart wiederanordnet, dass ein vorgegebener Progressionsmodus des Umfangs 

55 der Maschen erhalten wird, 

8. Verfahren nach einem der vorhergehenden Anspruche, dadurch gekennzeichnet, dass man zu einem emeuten 
Zusammensetzen aller Blocke ubergeht. 
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Claims 

1. Method for generating a 3D mesh in a structure whose geometry is known and imposed a priori, such as for 
example a geological structure defined from available data acquired by measurements, by analyses and/or inter- 
pretations of seismic data or logs produced in bore holes through the structure, with the object of producing a 
representative model of this structure, comprising prior modelling of the limits by triangulated surfaces, enveloping 
these limits as well as possible, and splitting of the zone into blocks having at most six faces, the faces of which 
correspond to previously triangulated major faults or discontinuities in the zone, the method being characterised 
In that it comprises: 

subdivision of each block into a set of elementary cells necessary to model the zone by finite elements, by 
creating a 3D lattice or mesh from each of the different blocks obtained, said subdivision comprising: 

a) creating a mesh from each of the faces of the block by interpolation from its edges so as to create an 
15 intermediate grid, followed by projection onto the triangulated surface of the nodes of this intermediate 

grid in a defined direction so as to generate face grids; 

b) population of each block with internal nodes by interpolation from the nodes of the different grids formed 
on the different triangulated faces of each block so as to obtain a uniform distribution of the nodes in the 

20 volume of each block; and 

c) rearrangement in the positioning of the nodes of the elementary cells created obtained by relaxation in 
the positioning of the nodes as a function of predefined criteria. 

25 2. Method according to claim 1 , characterised by successively a first operation for relaxation in the positioning of 
the nodes on each of the different triangulated faces, the population with internal nodes being effected from the 
relaxed grids, and then a second operation for relaxation in the positioning of these internal nodes. 

3. Method according to the preceding claim, characterised in that the first relaxation operation is carried out by 
30 iterative phases of relaxation and projection. 

4. Method according to one of claims 1 or 2 or 3, characterised in that a direction perpendicular to the intermediate 
grid is chosen as the direction for projection. 

35 5. Method according to one of the preceding claims, characterised in that a mode of interpolation of transfinite type 
is used. 

6. Method according to one of the preceding claims, characterised in that the positioning of the nodes of the ele- 
mentary cells created is rearranged so as to obtain uniform meshes. 



40 



45 



7. Method according to one of claims 1 to 5, characterised in that the positioning of the nodes of the elementary 
cells created is rearranged so as to obtain a given mode of progression for the size of the mesh. 

8. Method according to one of the preceding claims, characterised in thataW the blocks are reassembled. 
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Fig. 1 
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Fig. 2 
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Fig. 3 
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Fig. 6 
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